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Abstract 


Electromagnetics forms the foundation for electrical and computer engineering(ECE). Therefore, the first-principle solutions 
to ECE hardware-related problems and issues can be obtained by solving electromagnetic fields which are governed by 
Maxwell’s equations. Since only few analytical solutions are readily available for structures of regular shapes, numerical 
methods must be employed; many of them have been developed so far including the popular finite-different time-domain 
method (FDTD) and finite-element method (FEM). In the past 30 years, thanks to the tremendous increase of computer 
power, many software packages, based on different numerical methods, have been developed and made available in markets. 
Still, the challenges in computational memory and efficiency remain, in particular for Multiscale and Multiphysics modeling 
and simulations. 





In this presentation, we provide the summary of the three recent developments pertaining to the time-domain methods for 
electromagnetic modeling; they are the analytical FDTD solutions, a practical time-reversal method and a unifying node- 
based computational platform. The first development gives an in-depth look of the FDTD solution structures in terms of 
numerical eigen-mode expansions; it presents the ways to achieve multi-scale and multi-grid modeling. The second 
development allows the effective source locating and possible structure synthesis that avoids false solutions and can be 
carried out in a large computational domain. The third development shows the possibility of having most of numerical 
methods under a single mathematical framework and then hybridizing them. These three developments present potential new 
horizons for efficient and effective numerical modeling of Multi-scale and Multiphysics problems. Numerical examples are 
given to illustrate the points. 


Index terms: finite-difference time-domain (FDTD) method, finite-element method (FEM), time-domain modeling, eigen 
modes, multiscale, Multiphysics, time reversal, source reconstructions, node-based mesheless methods, unification, 
hybridization. 
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Introduction 


> Electromagnetic model: Maxwell’s 
equations of 26 equations by James 
Clark Maxwell in 1864; they quantify 
and summarize all the previous and 
experimental results and describe the 
behaviors of electromagnetic fields. 





James Clerk Maxwell. 


> Modern Maxwell’s Equations: 
simplification of the 26 equations to 
the 4 equations we see today by 
Oliver Heaviside in 1885 with 
introduction of flux and vector 
operations. 
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Introduction 


> Electromagnetic Fields: foundations of electrical and 
computer engineering > electric power, cars, radios, 
cellular phones, WiFi, computers etc. 


> Electrical and Computer Engineering Problems: solutions 
of Maxwell's equations: 


OD(x, y,Z,t) 
Ot 
OB(x, y,Z,t) 







= Vx H(x, y,z,1)-J(x, y, 2,1) 


= V x E(x, y, 2,1) 
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Introduction 


> Solutions of Maxwell's equations: 


e Exact analytical solutions: 
only possible for a few regular geometry of structures. 


% Approximate solutions: 


L Circuit theory and lump element models (approximated 
models of Maxwell’s equations): workable at low 


frequencies (i.e. structure size 10 times smaller than 
wavelength): Ohm’s Law, Kirchhoff’s Voltage & Current 
Laws. 
_} Numerical methods: 
e Time-domain methods; 
e Frequency-domain method. 
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Introduction 


> Time-domain methods: Three major techniques: 
% Finite-difference time-domain (FDTD) methods: 
o Derivatives in Maxwell's equations > finite-differences in Maxwell’s 
equations. 
% Finite-element time-domain (FEM-TD) methods: 
o Solution domains > discretized with triangular or tetrahedral 
elements; 
o solutions > expansions with linear shape functions in in each 
element; testing with Dirac impulse functions; 
o temporal derivatives > finite difference. 
“¢ Method of Moments - Time-domain integral equation (TDIE) methods: 
o Maxwell's equations > integral equations to be solved; 
o Solution domains > discretized with surface elements; 
o solutions > expansions with basis function in in each element; 
testing with pre-selected weighting functions; 
o temporal derivatives > finite difference in most cases. 


DALHOUSIE | 
Y UNIVERSITY Email: z.chen@dal.ca 








X1'an, China Slide 10- Oct. 17, 2017 


Introduction 


> Major Software packages: 
+ XFDTD: 


“CST 
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Analytical FDTD Solutions 


Eigen-mode decomposition of the FDTD formulation: 





> Discretization in space > 
Discrete Space 








= =-VxE 

at 

D _v.H_y Finite-difference [H]" 2-At[Dg][E"] 

Ot 1 
V.B=0 [E]” + At[Dy][H]""2 —= 
V:-D= oa 
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Analytical FDTD Solutions 


[HJ*= [H¿—At[Dg][E"] 


At 


[EJ += (EJ? + At[Dy] [H] — py" 





FDTD solutions: superposition of spatial 
eigen-modes modulated with time varying amplitudes 


Eigen MA vs 
spa MA eigenmode: | 
Stable mode : Vix 4 
Unstable mode: y; = 4 


Stable FDTD with large time step: discard unstable modes 
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Analytical FDTD Solutions 


Eigen-mode decomposition of the time-domain methods: 
take FDTD as an example 





FDTD solutions: superposition of spatial 
eigen-modes modulated with time varying amplitudes 





> To find eigenmodes [V|;, model-order reduction techniques 
maybe used since stable modes are of low-frequency 


e.g. run a few FDTD iterations 
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Analytical FDTD Solutions 


Eigen-mode decomposition of the time-domain methods: 
take FDTD as an example 










RTN 





Conventional 


The above 
FDTD applied | stable FDTD 
in coarse grid Applied in 
x fine grid 
E Ez | 8 
E MA 
i Only spatial 
interpolation is needed 
r-l [ 
y 2 E + Ey 
y y 





Stable FDTD solutions! 
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Analytical FDTD Solutions 


Eigen-mode decomposition of the time-domain methods: 


Sl > > > 


sr SUbgrid mesh becomes finer ==> 








M. Zhou, Z. Chen and W. Fan, “A subgridding 
scheme with the unconditionally stable explicit 
FDTD method,” in Proc. IEEE Numerical 
Electromagnetic Modeling and Optimization 
(NEMO), July, 2016. 
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Analytical FDTD Solutions 


Eigen-mode decomposition of the time-domain methods: 
take FDTD as an example 


1. Z. Chen and Peter P. Silvester, “Analytic solutions for the finite- difference time- 
domain and the transmission-line-matrix method”, Microwave Optical Tech. Lett., 
John Wiley &Sons, vol. 7, no. 1, pp. 5-7, Jan. 1994. 

2. M. Gaffar and D. Jiao, “An explicit and unconditionally stable FDTD method for 
electromagnetic analysis,’ [EEE Trans. Microwave Theo. Tech., vol. 62, no. 11, pp. 
2538-2550, Nov. 2014. 

3. Z. Chen, W. Fan and S. Yang, “Towards the wave-equation based explicit FDTD 
method without numerical instability”, in Proc. IEEE Intl Conf. Comp. 
Electromagnetics (ICCEM), pp. 265-267, Feb. 23-25, 2016. 

4. M. Zhou, Z. Chen and W. Fan, “A subgridding scheme with the unconditionally 
stable explicit FDTD method,” in Proc. IEEE Numerical Electromagnetic Modeling 
and Optimization (NEMO), July, 2016. 

5. J. Yan and D. Jiao, “A symmetric positive semi-definite FDTD subgridding 
algorithm for arbitrary grid ratios with uncompromised accuracy,” in Digests of 
IEEE International Microwave Symposium, June 4-9, 2017 Honolulu. 
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Analytical FDTD Solutions 


Analytic FDTD Solutions: 


FDTD solutions: superposition of spatial 
eigen-modes modulated with time varying amplitudes 





M sin(n—k)@ 

a” = > be (n—k 
2. x sino 

0. = arctan 





Source coefficient 
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Analytical FDTD Solutions 








Analytic FDTD Solutions: 


—+— FDTD method 
~*~ Analytic FDTD Solutions 


O 
A 


O 
N 


O 
O 





Electric Fleld Value (V/m) 
O 
N 


900 1000 1100 1200 1300 1400 
Iteration Number 


W. Fan, Z. Chen, and S. Yang, “On the analytical solution of the FDTD,” 
IEEE Transactions on Microwave Theory and Techniques, Vol.64, No. 11, Nov. 
2016, pp.3370-3379. 





SE DALHOUSIE 


UNIVERSITY Email: z.chen@dal.ca 





X1'an, China Slide 20- Oct. 17, 2017 


Outline 


Introduction 
Analytical FDTD Solution 


A Practical Time-Reversal Scheme 


A Unifying Computational Platform 


Conclusion 


DALHOUSIE | 
a UNIVERSITY Email: z.chen@dal.ca 





X1'an, China Slide 21- Oct. 17, 2017 








A Time Reversal Method 
Time-reversal and structure synthesis 


Normal forward wave propagation: 


Three probes recording 
the field signals 





~ Source location 
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A Time Reversal Method 
Time-reversal and structure synthesis: 








Reverse (backward) propagation: 


The probes take the recorded signals, 
time reverse them and reject them 
into the regions > backward simulations 
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A Time Reversal Method 


Time-reversal and structure synthesis: 








Simulation 


Ca ØST 


Forward Simulation Backward Simulation 
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A Time Reversal Method 


Time-reversal and structure synthesis: different source locations 











Failure 
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A Time Reversal Method 


Time-reversal and structure synthesis: the proposed formulation as an 
optimization problem - regularized least square (RLS) formulation 





N 2 


minimize ( N la|- Y y [7 | + Alal) 


\ n 
regularization 


Proposed Source Forward factor 
backward amplitude simulation 

simulation responses 

responses 
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A Time Reversal Method 


Time-reversal and structure synthesis: the proposed optimization 
formulation > results with impulse responses: full frequency band output 


signals 











or synthesized 
sources 
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A Time Reversal Method 


Time-reversal and structure synthesis: optimization formulation; in practice: 
band-limited output signals 




















Magnitude _ 
40} 
vi sal MN | dem Full-frequency-band: 
dad JUL dance. magnitude 
40! i i i 
es sel «a Limited-frequency-band: 
l nia dia , : | magnitude: practical signal 
O.O 0.5 1.0 1.5 2.0 2.5e11 
Frequency Hz 
40, Cbass Segre. | 
7 OPK PN A Pro, | Full-frequency-band: 
( | LAA a "wy e a, pp i | fm, qe lu | 
ao} ~“/ ww phase 
0.0 0.5 1.0 1.5 2.0 
40 


Limited-frequency-band: 
phase: practical signal 


cd) 





O.O 0.5 1.0 1.5 2.0 2.5 e11 
Frequency Hz 
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A Time Reversal Method 


Time-reversal and structure synthesis: optimization formulation; band- 
limited output signals; results. 
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A Time Reversal Method 
Time-reversal and structure synthesis: references. 


e W. Fan, Z. Chen, and W. Hoefer, “Source reconstruction from wideband 
and band-limited responses by FDTD time reversal and regularized least 
squares,” available online, IEEE Transactions on Microwave Theory and 
Techniques, August 2017 

e W. J. R. Hoefer, "Computational Time Reversal—A Frontier in 
Electromagnetic Structure Synthesis and Design’, JEEE Trans. 
Microwave Theory Tech., vol.63, no.1, Jan. 2015. 

e R. Sorrentino, L. Rosseli, and PB Mezzanote, “Time-reversal in finite 
difference time domain method”, IEEE Microw. Guided Wave Lett., vol. 
3—11, no. 11, pp. 402—404, Nov. 1993. 

e W. Fan and Z. Chen, "A condition for multiple source reconstructions 
with the time-reversal methods,” in Digest 2016 IEEE MTT-S 
International Microwave Symposium (IMS), San Francisco, CA, 2016, pp. 
1-4. 
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A Unifying Computational Platform 


A generalization theorv and a unifving computational platform: 





Node-based Meshless 











Circular 
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A Unifying Computational Platform 


A generalization theory and a unifying computational platform: 


The Method of Weighted Residuals (MWR) 
Node-based Meshless Method <4 Unifying Platform 












Method of 
Moment 
(MoM) 


Finite Element 
Method 
(FEM) 


Finite Difference 


Method 
(FDM) 
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A Unifying Computational Platform 


A generalization theory and a unifying computational platform: 





Node-based Meshless 


E(r,t) or H(r,t)= Y (rra, (t) 
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A Unifying Computational Platform 


A unifying platform: 
Node-based Meshless > Special case: Finite Element Method (FEM) 


C A 
E(r,t) or H(r,t) = YA r)a, (t) = 25 (rr, Ja (1) 


1.0 


NNs (r,r A x function 





Node positions are determined by triangular grids 
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A Unifying Computational Platform 


A unifying platform: 
Node-based Meshless > Special case: FEM > Special Case: FDTD 


p SR 
E(r,t) or H(r,t)= > b(n, )a,(0 = Tx), Oy YT, (27,94, 0 


p=! i.j.k, 


1.0 
JY 7 (r,r,) - Triangular or rooftop function 





Node positions are determined by rectangular grids 
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A Unifying Computational Platform 


A unifying platform: 
Node-based Meshless > Special case: FEM > Special Case: FDTD 


l. Z. Chen and M. Ney, “Method of weighted residuals: a general approach to 
deriving time- and frequency-domain numerical methods,” /EEE Antennas and 
Propagation Magazines, February, 2009 (a feature article), pp.51-70 

2. Z. Chen and S. Luo, “Generalization of the finite-difference based time-domain 
methods using the method of moments,” IEEE Trans. on Antennas and 
Propagation, Vol. 54, No. 9, September, 2006, pp. 2515-2524 

3. Z. Chen and Y. Yu, “Node-based methods: a unifying approach to developing 
numerical techniques,” invited special session paper, in Proc. of Intl. Conf. on 
Electromagnetics in Advanced Applications, Sept. 12-16, 2011, Torino, pp.737-740 

4. J. Wang, Z. Chen, et. al., “A computational platform that unifies and hybridizes 
numerical methods,” in Digest 2017 2017 IEEE International Symposium on 
Antennas and Propagation, July 9-14, 2017, San Diego 
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A Unifying Computational Platform 


A very good but unfortunately not-well known early manograph by 
Prof. Wen-xun Zhang implies some of the ideas of the unification of 
the different numerical methods 


"ENGINEERING 
ELECTROMAGNETISM: 
FUNCTIONAL 

METHODS 


y 
Wen-Xun Zhang É 





W. X. Zhang, Engineering Electromagnetism: Functional Method, 
Ellis Norwood, 1991 
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A Unifying Computational Platform 


Multiphysics and multiscale modeling 
Example: electric fields equation interact with the thermal equation? 


E° 
Electric fields equation ==> P(x, y) Source of the heat 


J 


Heat transfer equations I 








Parameters will be changed with 7 
l 


Pid +aT (x, y)) 





o(x, y) = 
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A Unifying Computational Platform 


Multiphysics and multiscale modeling 
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A Unifying Computational Platform 


Multiphysics and multiscale modeling 


Simulation 





Dirichlet boundary condition Robin boundary condition 
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Conclusion 


Eigen-mode and Analytic Representation of 
Numerical Solutions 


> The FDTD solutions are found to be the superposition of the 
Spatial eigen modes. 


> How to deal with these spatial eigenmodes may hold the key 
for further improving the FDTD computational efficiency. 


> Analytic FDTD solutions have been found, promising for 
efficient time-domain modeling. 


> Similar can be applied to other time-domain methods. 


>» Non-conformal or unstructured methods are naturally 
effective for multiscale modeling. 
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Conclusion 
The Proposed Time Reversal Method 


> Time reversal techniques may be an alternative way 
to synthesize the sources or structures if the forward 
simulation results can be specified as the desired 
responses. 


> A optimization formulation with the time reversal 
method has been proposed as an effective way to 
reconstruct sources without spurious results. 


>The proposed method works in the case of 
bandlimited signals at the output nodes, making the 
method practical. 
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Conclusion 
A Unifying Computational Platform 


> There is not a single numerical method that is 
good at solving all the electromagnetic problems 
(e.g. multiscale or Multiphysics) interfacing or 
hybridization of different methods may solve them 
all > a unifying platform is then desirable. 


» The node-based mesheless method can be used 
as the unifving platform although its theory and 
applications have not been fully established vet. 
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Congratulations and Happy Birth 
to Prof. Wen-Xun Zhang! 
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